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Examples of glasses are abundant, yet it remains one of the phases of matter whose understanding 
is very elusive . In recent years, remarkable experiments have been performed on the dynamical 
aspects of glasses. Electron glasses offer a particularly good example of the 'trademarks' of glassy 
behavior, such as aging and slow relaxations. In this work we review the experimental literature on 
electron glasses, as well as the local mean-field theoretical framework put forward in recent years to 
understand some of these results. We also present novel theoretical results explaining the periodic 
aging experiment. 
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I. INTRODUCTION 

Slow relaxations in nature have been observed long ago, in many contexts. Two important examples are the 
'stretched exponential' behavior experimentally observed by Kohlrausch in the 19th century, when studying discharge 
in a certain electrostatic setup known as a Leiden jar [T] and mechanical relaxations in silk threads, used by Weber 
[2] to hang his magnets in his seminal works on magnetism. 

Slow relaxations are also important signatures of glassy behavior: the Vogel-Fulcher law [5] states that the timescale 
for the thermalization of a glass diverges at a finite temperature. When valid, this temperature is the phase-transition 
to the glassy phase, below which non-ergodic behavior occurs: the system does not equilibrate, and some states are 
never visited. Slow relaxations have been known to occur in magnetic materials known as spin glasses @HI]> vortices 
in superconductors in the glassy phase [8], and in the electron glass system, to name but a few. 

In many cases in nature the slow relaxations are logarithmic, over several decades in time. These systems range 
from flux relaxation in superconductors [S] , through crumpling paper |10j where the volume is measured as a function 
of time, as well as biological systems such as plant rheologypT], and electron glasses which we will focus on here. 
One should note that it is not trivial to distinguish a logarithm from a stretched exponential: one has to go to 
short time scales, where the power-law behavior of the stretched exponential is distinctly different from logarithmic. 
For this reason relaxations in the electron glass were initially erroneously reported as stretched exponential. Fig. [T] 
demonstrates a logarithmic relaxation after a sudden change in the voltage of a gate coupled to the sample. 




t (sec) 



FIG. 1: Demonstration of a logarithmic relaxation in electron glasses, after a sudden change in gate voltage. In some cases, 
the logarithmic change in conductance can be measured from times of order of seconds to several days [T2] . Data courtesy of 
Z. Ovadyahu. 
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In this review we shall explain how one can understand these slow, logarithmic relaxations in electron glasses via 
the use of a local mean-field approximation, and also extend this to various aging protocols, in which the response of 
the system to a perturbation depends also on the time the perturbation was applied for. We describe the experimental 
details, but only briefly discuss theoretical frameworks other than the local mean-field approach. 

The structure of the manuscript is as follows. We first explain different mechanism leading to a broad distribution of 
relaxation times, which will give rise to logarithmic relaxation over a broad range. Wc review the different experimental 
protocols demonstrating aging, and their results, focusing on granular aluminum and InO. We then go on to describe 
the local mean-field approach, and compare its prediction with the experimental results. We discuss the possible 
relation between 1// noise and logarithmic relaxations. Next, more involved protocols related to memory effects in 
the system are reviewed, and the necessary theoretical framework for a qualitative understanding of them is given. 
We conclude with our view regarding the important, open questions in this field. 



One may decompose the logarithmic relaxation into a weighted sum of decaying exponentials, by taking an inverse 
Laplace transform. This gives a distribution of relaxation rate which is P(X) = c/A, with c = l/log(A max /A m i n ), 
with A max and A m i n the upper and lower cutoffs of the distribution. Notice that the distribution of the logarithm 
of A is uniform between the lower and upper cutoffs. Indeed, taking its Laplace transform of P(A) yields, for times 
1/A max < t < 1/A min Q3]: 



with the Euler constant. Notice that the upper cutoff does not enter the expression, but the lower one does. 

In order to understand the logarithmic relaxations, it is sufficient to understand how a distribution P(X) = C/A 
arises. This is the purpose of the next section. 



The simplest way to obtain a broad I/A distribution, is via thermally activated processes: let us assume that a 
particle hops between trapped states, and that the distribution of the energy barriers between two adjacent traps is 



II. WHY A LOGARITHM? 




(1) 



A. Activated processes 



uniformly distributed in the interval [E m [ n , E max ] , P(E) 



l 



If we assume thermally activated processes, 



E, 
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the rate A associated with a given barrier is given by the Arrhenius law, A ~ er E l T (throughout the manuscript 
Ub = 1)- Calculating the distribution of the rates, we find that P(X)d\ = P(E)dE, and therefore P(X) ~ T/X. 

In various experiments |14H16j the temperature dependence was shown to be much weaker, if at all measurable. 
For these cases, this rules out the possibility of the above mechanism as being responsible for the slow relaxations, 
and opts for a quantum mechanism, nearly independent of temperature. Quantum tunneling is the natural process 
to look at, in order to avoid the strong dependence on temperature. 

B. Tunneling processes 

If we look at a particle hopping between localized states, the tunneling rate is exponential in the distance r (up 
to polynomial corrections). If we assume that the spatial distribution of localized states is uniform, we can readily 
calculate the distribution of hopping rates between a given site and its nearest- neighbor: the rate A ~ e~ r /£, where £ 
is the localization length, taken here as a constant for simplicity. We have P(X)dX = P(r)dr, with P{r) the probability 
distribution of having the nearest-neighbor at a distance r. It is easier to calculate the cumulative of this distribution, 
i.e., the probability C(r) that all neighbors are more than a distance r away from a given site. The probability of 
a given site to be more than a distance r away is p = 1 — Vdr d /L d , with Vd — n d ^ 2 /T(d/2 + 1) the volume of a d 
dimensional unit sphere and L the system size. Thus, C(r) = p , where N is the number of sites. In the limit N 3> 1, 
this can be simplified to give: 

P(r) = ^[r/(r)] d -V^/<< (2) 

( r ) = pfi/d being the average nearest-neighbor distance. 
This immediately gives: 

P( A ) = d MM [-flog(A)/(r)] d -V v *'-* to gW/< r >r (3) 
A 

Thus, we see that this mechanism yields a 1/A distribution, up to logarithmic corrections, which may be of impor- 
tance if the relaxation is probed over enough decades in time. Related analysis leading to approximately logarithmic 
relaxations was made in Refs. [17Ul9j . where the relaxation rate was assumed to be the exponential of a smoothly 
distributed variable. 

Ref. |20j analyzes this problem more carefully, and finds the exact distribution of relaxation rates, for any spatial 
dimensions, in the low density limit. It is shown that the resulting distribution is approximately P(X) ~ 1/A, but in 



dimensions higher than one there are logarithmic corrections similar to those of Eq. Q , but with different numerical 
coefficients. The real-space renormalization group method used also allows to find the structure of the eigenmodes, 
which show interesting localization properties: all eigenmodes are localized, but with size diverging as one goes to 
vanishing eigenvalues, as ~ e c ' log (~ A / 2 )I, with C a constant. In one dimension, for example, this would yield a 
power-law relation. 

An interesting and relevant question addresses the issue of how this picture is changed when one considers tunneling 
not only of a single particle, but simultaneous quantum tunneling of a number of electrons 21 -28 . Obviously, the 
associated timescales would be longer, since now one has to replace e~ r ^ by e — r ^'£, where Tj are the tunneling 
distances of each of the particles involved in the many-particle tunneling process. In principle there are many pos- 
sibilities to connect the initial and final states, but the one where the sum over distances is minimal would be the 
dominant one. There should be quantitative changes in the form of P (A), namely, in the corrections to the numerator 
of Eq. Since estimating the numbers for single-particle processes (see section (VIE)) yields too short relaxation 



times, we are led to believe that these processes must be considered. A step in that direction was taken in Ref. [29) . 
However, the temperature dependence which they estimate is not consistent with the rather temperature independent 
results discussed above. In our view, the question of many-particle tunneling still deserves further attention. 



III. AGING EXPERIMENTS 

The first experiments showing slow relaxations in the electron glass were performed by Monroe et al. (3U] in 1987, 
who looked at the time-dependence of the capacitance of a sample of GaAs, after injection of excess charge. 

Later, Ovadyahu and collaborators have shown remarkable aging behavior (a concept to be shortly explained) in a 
series of extensive experiments done mainly on InO samples, in a field effect transistor setup [HI HH [T51 Q~71 Qj|l |3T] . In 
the following, we review the basic experimental procedure used in the experiments, and clarify a common confusion 
between two different kinds of aging experiments used by the scientific community We then go on to review the 
existing theoretical models used to explain the results. 



A. Aging protocols 



Aging is the general phenomenon related to a relaxation which depends on the 'age' of the system: in some 
experiments, the system is quenched to a low temperature phase at time t = 0, and the response R of the system 



to a perturbation applied at time t w is tested at time t. If R depends explicitly on both t and t w , and is not only a 
function of t — t w , the system is referred to as "aging" over time, since its properties are not translationally invariant 
with respect to time. 

Indeed, experiments done on spin glasses [3HZ] have shown such "aging" behavior. In this context, the above 
protocol is known as a thermoremanent magnetization (TRM) experiment. 

Another protocol, which we will focus on, follows the following procedure, illustrated in Fig. [2} 




t (sec) 



FIG. 2: A schematic description of the different stages of the IRM aging protocol, taken from Ref. [32] .The lower part shows 
the experimental result (circles), courtesy of Z. Ovadyahu, and theoretical predictions of the model (solid line), discussed in 
section fiV) . Stage II of the experiment shows exactly the same experimental data as Fig. [TJ i.e., at this stage the conductance 
relaxes logarithmically with time. In stage III, the relaxation turns out to depend only on the ratio t/t w . 



1) The system is let to equilibrate for a long time t\ at the low temperature phase (not necessarily reaching the 
true equilibrium). 

2) A perturbation is applied for a time t w . 
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3) A time t after the perturbation has been switched off, the physical observable f(t,t w ) is measured. 

It is important to note that in practice, in most cases the system does not fully equilibrate during the first stage, 
and in that case this timescale can have implications on the observables in later stages. 

The main experimental result for electron glasses are as follows: During stage II, the relaxation is well described 
by a logarithm, while in stage III excellent data collapse is obtained when time is rescaled according to t w , in other 
words, the relaxation f(t,t w ) in fact only depends on the ratio t/t w . 

In the spin glass context, this type of experiment is known as an isothermal remanent magnetization (IRM) ex- 
periment. Fig. [3] illustrates schematically the two protocols described above (TRM and IRM), as well as two other 
protocols which will be discussed later on. 

In a sense, the IRM experiment is much simpler, since the system is presumably close to equilibrium at all times 
during the experiment, provided the perturbation is small such that we are in the linear-response regime. Recently, it 
has been suggested that calling this type of experiment an aging experiment is misleading, exactly for this reason |33] 
(since this experiment does not show that the system properties depend on the time from the thermal quench). As 
we shall later see, understanding the slow relaxations also for this experiment is actually quite involved, and presents 
some unsolved questions. 

We shall now review the experimental results obtained for the different systems. 

1. Indium Oxide 

Pioneering experiments using the above IRM protocol, in the context of electron glasses, have been performed 
mainly on Indium Oxide, by Ovadyhau et al. In Ref. [14] they show that "full" (also called simple) aging is obtained: 
the function f(t,t w ) depends only on the ratio t/t w , i.e., f(t,t w ) = g{t/t w ). For short times, they observe that 
the function g{t/t w ) depends on the logarithm of its argument. These properties will be elucidated in the next 
subsection. It turns out that this is by no means an accident related to the peculiar properties of InO, and occurs also 
for various other systems. Qualitatively similar results are obtained for crystalline samples and amorphous samples, 
demonstrating the broad applicability of the results. Ref. [TU] discusses the deviations of f(x) from a logarithmic 
dependence, and also observe deviations from full aging for large enough gate voltages. In Ref. [12], Orlyanchik et al. 
use the same protocol but perturb the system by an electric field (which they termed 'stress aging'), and not a gate 
voltage. This did not modify the full aging behavior, for not too large fields. 

Fig. [4] demonstrates this scaling behavior for a number of cases discussed in the following, by showing data collapse 
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FIG. 3: Schematic illustration of different protocols in use, relating to the aging and memory effects in electron glasses. In 
all protocols, at a given time (the crossing point of the x and y axis, in the figure), the temperature is quenched to reach the 
glass phase, (a) describes the TRM protocol, only recently performed in the context of electron glasses [33], but extensively 
used in spin glass experiments [JH?]- In this protocol, one waits a time t\, and perturbs the system (in the electron glass, 
typically by changing the gate voltage). A time t later, the response of the system (i.e., the conductance) is measured. If 
ii 2> t, one obtains a logarithmic relaxation of the conductance, for times t smaller than the cutoff 1/A m i n of the relaxation 
rate distribution. Such a relaxation is shown in Fig. [I] The important point, however, is that generally one can also consider 
ti < t, which distinguishes this from the IRM protocol. The measured conductance will depend explicitly on both variables ii 
and t, and in a non-trivial way, which is not well understood, (b) describes the IRM protocol, in which t\ is assumed much 
larger than other experimental timescales involved (in typical experiments, it is of the order of a day). After the time t\ one 
changes the gate voltage, for a duration t w (the 'waiting-time'). A time t after the gate voltage is returned to its initial value, 
the conductance is measured. Here, the ratio t/t m can be much smaller or much larger than unity, as long as both are still small 
compared to t\. Fig. [2]shows the experimental results of a typical experiment, and a curve showing the theoretical prediction. 



The theoretical analysis for this protocol is performed in section ( IV A 2 1 . (c) describes a generalization of this, where a periodic 
square pulse gate voltage is applied, and the conductance is measured throughout the experiment. Unlike most systems, here 
the response to a periodic signal is not periodic, since the whole experiment is essentially still in the transient period, before the 
periodic regime is maintained. Fig. [5]shows the result of such an experiment. The relevant theoretical analysis is performed in 



section (IV A 3 1. It is assumed that t <C t\. (d) describes the two-dip protocol: here, the gate voltage is changed at some time, 



and the system begins to relax logarithmically to its equilibrium. A time t later, a scan of conductance versus gate voltage is 
made. It turns out that as the system equilibrates at some gate voltage V , it 'digs' a dip in conductance at this gate voltage. 
In this protocol a new dip begins to form at the new gate voltage, as the old one is erased over time, which gives the protocol 
its name. The protocol is used to demonstrate the memory effects, and also to quantify the timescales involved. It is described 



in section I VB|, and analyzed qualitatively from the theoretical point of view in section (VID|. Also here it is assumed that 
t < h. 
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when the time axis is rescaled by t w . The theoretical curve will be discussed in section ( IV I . 



2. Granular Aluminum 



Extensive experimental work using the above IRM protocol has been performed on granular Aluminum by Grenet 
et al. [El [34]. Most of the results show striking similarity to those obtained for InO, see Fig. [4] Some differences 



occur in the behavior of the two-dip experiment, which will be described in section (VB|. 



3. Other materials 

Experiments performed on a 2D electron system in silicon by Jaroszyhski and Popovic have shown slow relaxations 
and aging |35H37) . They experimentally observe a transition between a phase where aging effects exist (dependence 
on t w ) to one at higher densities where no such dependence exists. In the aging regime, they find full aging only 
below another critical density. While very interesting by its own right, it seems that the underlying mechanism is 




FIG. 4: Comparison between the theoretical prediction of Eq. (17 1 and various experiments. No fitting parameters are used. 
Data courtesy of Z. Ovadyahu and T. Grenet. 
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different than the one which occurs for InO and Al, which is the one we focus on in this review, and therefore we shall 
not elaborate on this extensive set of experiments. Borini et al. |38) did observe full aging in porous silicon, at room 
temperature, which seems to fit the theoretical framework that we shall put forward in section ( IV ) , see Fig. [4] 



In the group of A. Goldman, experiments performed on thin films of bismuth or lead showed logarithmic relaxations 
in the conductance after a change in gate-voltage, which relaxed to a different conductance value, depending on the 
value of the gate voltage [39] ■ They also observed an anomalous field effect [3D], which will be discussed in section 
Q. It would be interesting to see the behavior of this system when the IRM protocol is applied. 

Recent experiments have shown that Nickel samples can also exhibit slow relaxations, with behavior which seems 
to be similar to that of granular Aluminum and InO (41) . In this case, however, a transition between a ferromagnetic 
state and a super paramagnetic state seems to have a substantial effect on the relaxation rates. As a result, the 
magnetic held has a striking effect on the relaxation times. Future experimental and theoretical work in this direction 
looks promising. 

B. TRM in electron glasses 

In a recent experiment Grenet et al. have extended the above protocol to the case when the system is only partially 
equilibrated before the measurement sequence is performed [33] . This adds another timescale to the problem t\ 1 the 
time from the temperature quench to the first voltage step, as shown in Fig. [3]ja). The above IRM protocol is the 
case where t\ ^> t w . It turns out that for smaller t±, the measurements depend on it. Remarkably, a 'superposition 
principle' akin to that observed in spin glasses |42j is still maintained: the response to two voltage pulses is the sum 
of the responses to each of the them independently. This has not been explained theoretically, for the spin glass nor 
the electron glass, and would be an interesting subject for future research. 

C. Generalized IRM (periodic protocol) 

Another protocol that comes to mind is the generalization of the IRM protocol to deal with an arbitrary series of 
pulses. An elegant experiment is discussed in Ref. |43j . where the sample has been subjected to a periodic sequence 
of square pulses, see Fig. [3jc) . Unlike the common scenario, where the response to a periodic signal is also periodic, 
after a short transient, here, the whole experiment is found in the transient period. Fig . [5] shows the result of the 



experiment, together with the theoretical prediction, which will be derived and explained in section (IV A3) 
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IV. WHY IS FULL AGING OBSERVED AND WHEN DOES IT FAIL? 

It is the purpose of this section to present the theoretical models explaining the aging behavior, and the limitations 
to obtaining full aging. 



A. Local mean field approach 

In a pioneering work, Thouless, Anderson and Palmer introduced the TAP equations, to deal with the equilibrium 
properties of spin glass [44]. A similar approach was used in Ref. [45], to deal with the statics of the Coulomb 
glass problem. It was shown that by using a local mean-field approach (which retains the individual identity of each 
site - not averaging over the disorder), one obtains the Coulomb gap, a soft gap in the density-of-states discussed in 



section ( VI A 3 1 . The long range nature of the problem gives rise to a large number of effective neighbors a given 
electron interacts with, which gives intuitive justification for the use of the mean-field approach. This was quantified 
in various works |46H48] . indicating that the relevant criterion for the validity of the mean- field is that the power a 
characterizing decay of the interaction l/r a should be small enough. The results suggest that for Coulomb interactions 
in two-dimensions and above one could use the local mean-field approach. 

In Ref. [33] j the local mean-field approach was extended to treat the dynamical aspects. 




Time (seconds) 



FIG. 5: Comparison between the theoretical analysis of section (IV A3 1 (solid line) and an experiment (circles) where the 
sample is subjected to a periodic perturbation. The dashed line demonstrates the logarithmic increase of the baseline, see Eq. 
No fitting parameters are used. Data courtesy of Z. Ovadyahu. 
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The equations describing the dynamics are: 



drii 
~dt 



(iji ~ Hi) , 



(4) 



with the rates given by Fermi's golden rule [50]: 



7ij ~ \M q \ 2 vfi(l - fj)e—¥-[l + N(AE) 



(5) 



where f{ is the Fermi-Dirac distribution, AE is the energy difference before and after the tunneling event, v is the 
phonon density-of-state, M q is the corresponding matrix element and N(AE) is the Bose-Einstein function. For 
upward transitions (Ej > Ei) the square brackets are replaced by N(AE). These rates may be renormalized due to 
polaron-type orthogonality effects [51] , which were suggested to be of importance in electron glasses [52] ■ 

The Bose-Einstcin function expresses the underlying assumption that the phonons thermalize fast and are found in 
thermal equilibrium. AE is the energy difference between the two sites, which is affected by the Coulomb interactions 
with all other sites, i.e., AE = Ej — Ei, with: 



It should be emphasized that in this local mean-field approach the sites are not assumed to be equivalent, as is 
clearly seen by the implicit site indices. This is very different than other mean-field approaches to the electron glass 
problem [53H56] , 

A characteristic of glasses is the abundance of metastable states, with energies close to that of the true ground state 
[57j . Also in this case, there are many metastable states, i.e., electronic configurations for which all the time derivatives 
(the LHS of Eq. Q) vanish. This has not been investigated systematically within this framework, although it would 
be interesting to do so. A numerical investigation of the energy landscape has been performed in a related model 
through Monte-Carlo simulations, by Baranovskii et at, where the metastable states were coined "pseudo ground 
states" [58] . For recent related works see [59l |60] 

Close enough to a particular metastable state characterized by an electron configuration no, one can linearize the 
equations of motion for 5n = n — uq: 




(6) 



(7) 
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where the off-diagonal matrix A is given by the expression 



^-^•„0(i_ n 0) E 



1 , 1 1 



(—-—). ( g ) 



with 7y the equilibrium rates obeying detailed balance, see Eq.(j5J). 

The diagonal matrix elements are given by An — — A^j , guaranteeing particle number conservation. The 



matrix is real but not hermitian, due to the n>j(l — rij) term. In section (IVBll we shall give a very natural 



explanation for the explicit form of A, related to the properties of the noise via the Onsager relations [61]: 



A = 7 /?, (9) 

where 7 is a symmetric matrix describing the transition rates at equilibrium, sec Eq. ([5|, and with the diagonal 
defined such that the sum of columns vanishes, and /3 is a symmetric matrix whose inverse gives the equilibrium 
(equal-time) correlation function between two given sites: 

^ = ^W) + r 2 i^- (10) 

While A is not symmetric, its eigenvalues are nevertheless real: this can be seen using Eq. Q 6 1]. Near a stable 
minimum, the eigenvalues must be be negative |49j . and their distribution will determine the relaxation properties. 

1, Properties of the relaxation matrix A 



In section ( |IIB ), the statistics of nearest-neighbor distances was analyzed, and was shown to be distributed as 



P(A) ~ 1/A up to logarithmic corrections, due to the exponential nature of the tunneling processes. Looking at 
the matrix A of the previous section, we find that the off-diagonal matrix elements also decay exponentially in the 
distance. A numerical study made in Ref. |49j showed that in spite of the more involved form of the matrix elements 
(related to the energy dependence), the above property is still retained. Considering the toy-model where this energy 
dependence of the matrix elements is neglected (which is physically correct when the localization length is small 
enough, but is not always the case in experiments), one obtains a random-matrix class which is defined as follows: 

1. Aij = e"W€. 

2. A i2 = — J2j^i Aij. 
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As mentioned earlier, the second property arising directly from particle conservation number. 

A heuristic approach to understand qualitatively the emerging distribution of relaxation rates was made in Ref . |49j . 
Later, an exact solution was found, in the low-density case |20j . which confirms the numerical result described above: 
up to logarithmic corrections which depend on the dimensionality and may be of importance if the experimental 
resolution is fine enough, the leading order behavior is a uniform distribution of the logarithm of the relaxation rate, 
i.e., P(X) ~ 1/A. This will play a key role in the understanding of the IRM protocol, to be described next. 



2. IRM protocol - theory 

In Ref. [32], the -P(A) ~ 1/A distribution was the starting point of an analysis of the IRM protocol. An underlying 
assumption is that the state of the system can be described by a vector Sn, describing the deviation from an equilibrium 
or metastable state, and that the evolution in time is given by: 



— = ASn. (11) 

Notice that although the above equation is identical to Eq. Q, which was derived from the mean-field approach, 
here we take a more general approach, where this equation is the starting point. The distribution mentioned above is 
that of the eigenvalues of the matrix. 

The crux of the matter is that when a gate- voltage, for example, is changed, the equilibrium point changes, say from 
ii a to fib, such that the system is instantaneously thrown out of equilibrium. Relaxation to the new equilibrium is 
excitation with respect to the initial equilibrium. Fig. [2] describes the different stages of the experiment schematically. 

It is clear that in the stage II, the system is relaxing back to its new equilibrium, so that in terms of the eigenmodes 
its deviation from equilibrium can be written as: 



n(t) = n b + J2 c q b q e- xb * ( t +t ~\-t va < t < 0, (12) 

where A* are the relaxation rates of the relaxation matrix A at point B, b q are the eigenvectors, and c q are the weights 
of the decomposition into eigenmodes. 

Let us assume that the modes contribute, on average, positively and uniformly to the conductance a(t), i.e.: 
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So ~ \ A > 



(13) 



where \A\\ is the amplitude of the mode A in the decomposition of ri(t) — tit, into eigenmodes. This is a central point 



which will be discussed in section (VI B). We obtain that 



Sa(t)~J2e- X « {t+t ™\-t w <t<0. 
i 

For the large time window 1/A max <t< 1/A m i n , this gives a logarithmic relaxation: 



(14) 



5a(t) ~ -7b - log[fA mi „], 



(15) 



with 7£ is the Euler constant, and t = t w + 1 is the time from the first change in gate voltage, obviously independent 
of t w . 

Similarly, in stage III, the relaxation is described by: 



= n a + y^c g a g ( 



e i tu 



l)e 



(16) 



An equivalent sum was given heuristically by Grenet et al. |16j . 

For the IRM protocol, we finally obtain after summing up the series in the same manner as that done to obtain 



Eq. (15) [32]: 



8a(t, t w ) oc log(t + t w ) - log(t) = log(l + t w /t). 



(17) 



Fig. |4] shows the correspondence of this to the experimental results in 3 different electron glasses. 



We are now in a position to understand why full aging occurs within this model: in general, Eq. ( 16 ) predicts a 



functional dependence on the parameters t and t w of the form g(t + t w ) — g(t). For the specific underlying distribution 



of relaxation rate, it was shown that g{x) is logarithmic, which leads to the full aging of Eq. (17). It can be shown 



that the converse is also true, and that P(A) ~ 1/A (associated with a logarithmic g(x)) is the only distribution which 
will give rise to full aging. 
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It is important to note that in the above analysis it was assumed that the system was initially in equilibrium. As 
mentioned, in practice it is possible that the system did not fully equilibrate in the first stage of the experiment, and 
in that case the reciprocal equilibration time 1/ii (see Fig. [3]) will replace 1/A min . 

3. Generalized IRM (periodic protocol) - theory 

It is possible to extend the IRM protocol, to the case where a more complex sequence of pulses is used. A natural 
choice is a periodic sequence of pulses, i.e., the system state is switched from one value of gate voltage to another, 



after a t w length of time. Such a protocol was carried out in detail in Ref. [S3] , as discussed in section ( III C I . We 
shall now give a theoretical analysis of the experiment, in the same fashion as was done in the previous section: 

At every instance in time, the system relaxes to its current equilibrium, which we shall denote by A or B, depending 
on the current value of the gate voltage. Every time the gate voltage is switched from A to B, the vector dn gains an 
extra An. Similarly, when the voltage is switched from B to A, the gain is — An. 

Thus, for example, a time t after two cycles of gate voltage changes, the state of the system is described by: 

n(t) =fi a +J2 c g a q [l - e A ^ (e"^*- - lJe-^-Je"^*. (18) 

Assuming, as before, that points A and B are close enough such that we can neglect the difference in the eigenmodes 
and eigenvalues, the sum can be readily evaluated to give: 



5<j ~ log[l + *„,/(* + 2t w )] + log[l + t w /t}. (19) 

Writing out the sum explicitly for a general sequence of pulses, we find the following simple rule: the contribution 
of a step up is log[A m i n (i — t ste p)], and the contribution of step down is — log[A m i n (t — t step )]. Indeed, for the IRM 
protocol we obtain log(A m i„(t m + t) - log(A min (i) = log(l + t w /t). 

Looking at the system state a time t after N pulses, its conductance is given by: 



5a ~ J2 N (-±y +1 log[* + *»(2j - 1)] - log[t + t w (2j - 2)]. (20) 

3=0 

To sum up the sequence, it is helpful to sum subsequent pairs of switches. Looking at the contribution of a pair 
of switches occurring a time t ^> t w earlier, it contributes log(l + t w /t) ~ t w /t. Therefore the response after N ^> 1 
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pulses is given approximately by: 



JV 

t w /(2jt w ) = log(iV)/2 = log[V(2^)]/2. (21) 

3=1 

The fact that this series diverges is profound: it means that we will have approximately periodic signal, but with 
a baseline that diverges logarithmically in time. This is exactly what is demonstrated in Fig. [5j comparing the 



evaluation of Eq. ( 20 1 , with no fitting parameters, with experiments . 

The striking result of this analysis, confirmed by the experiment, is that while we are used to systems responding 
periodically to a periodic perturbation after a short transient (e.g: RCL circuits, mechanical systems etc.), here, the 
system remains in the 'transient' period for the whole duration of the experiment. 



B. Noise in electron glasses 

1/f noise occurs in many physical [521 as well as biological and economic systems [5H[55]. Electron glasses 
are no exception, and various experimental studies have measured 1/f noise in such systems 66 - 72J . An even larger 
amount of work has been invested in the theoretical aspects of this problem. The pioneering works of Shklovskii |73j 
were later followed by various other approaches [2H1 ITiHTS] . In a recent work [5T] , the mean- field framework has been 
shown to yield a 1/f noise in the site occupancies, which has not been measured until now, but being a local object 
is a much simpler property than the conductance noise. 

In the following, we shall shortly explain the relation between the slow relaxations extensively discussed earlier, and 
the 1/f noise. We note that the correspondence to the experimental data should be done with extra care: since 1/f 
is so common, the experimental observation of such noise by itself is by no means a confirmation that the electron 
glass mechanism yielding this noise is the dominating one, and additional signatures must be considered. 



1. Why 1/f noise occurs? 

Although noise is measured at equilibrium, Onsager understood in a seminal work that it is generically related to 
the way a system returns to equilibrium after a slight perturbation |80) . 

In the following, we shall show how this Onsager principle relates the logarithmic relaxations discussed earlier to 
the noise in the system. 
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Onsager's regression hypothesis [81] states that the equation of motion of the correlation function <j>ij(t) 



(5ni(t)5rij(t)) is obtained by simply replacing 5rii(t) in the equation of motion (11 1 by the function (j>ij(t): 



dcj>jj{t) 
dt 



Aik4>kj{t)- 



(22) 



To find the correlation function one also needs the initial conditions: 



c/> ij (0) = (Sn i (0)Sn j (0))=pr j 1 . 



(23) 



To find the matrix j3, we write the free energy |45j : 



{\ + nt) log(- + hi) + {\- n%) log(- - hi) 



(24) 



with hi = ni — |. The local mean-field equations can be obtained from the minimalization condition = 0. Notice 
that each site contains a positive background charge of |, to keep charge neutrality. 
Expanding F near a metastable state (a local minimum), we have: 



kT 



F = F + — ^PijSruSnj, 



(25) 



with /3 given by Eq. ( 10 1 



The correlation matrix is proportional to /3 , since we have a quadratic free energy: 



{5nMj) = {p- l )i. 



(26) 



Regardless of the question of determining the noise spectrum, which we shall shortly derive, at this stage we can 
use another theorem due to Onsager to put the form of the non-hermitian relaxation matrix A of Eq. ([8| in a more 
natural context: the Onsager symmetry principle states that A(3~ l must be a symmetric matrix. Indeed, in our case 



A = 7/3, with jij the equilibrium current between sites i and j, as was mentioned in section (IV A) 



Let us proceed to the calculation of the noise spectrum. After finding /3, we can solve Eq. (22 1, and obtain the 
noise spectrum of the average site occupancy |61j . This leads to the expression: 



(27) 
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where ipf describes the amplitude of the a'th eigenmode of the relaxation matrix A. The factor J- describes the 
physically clear fact that only sites with energies close to the Fermi-energy will contribute to the noise, since they are 
"soft" and their occupancy fluctuates significantly. We therefore find that: 



a . t 



where ((; )) denotes averaging over sites as well as time. 

Plugging in the -P(A) ~ 1/A distribution of rates, we finally obtain the noise spectrum: 

{{Sn2)) " ~ h £7 dX ^w = ^ XT dm rh- (29) 

This shows that for A m j n <C oj <C A max , a 1/f spectrum indeed follows for the noise in the average occupation 
number. 

V. MEMORY IN ELECTRON GLASSES 
A. Anomalous field effect 

There is a clear experimental demonstration that for electron glasses perturbing the system causes the conductance 
to increase. In 'ordinary' field effect transistors, where interactions are not as strong, adding charge carriers to the 
system by changing a gate voltage (see Fig. [6]), will increase the conductance, and taking charge away will make it 
decrease. In electron glasses, however, after full or partial equilibration, changing the charge carriers number in any 
direction will make the conductance increase. This surprising property is called the anomalous field effect |82j . It 
is clear that the only thing which distinguishes the particular value of gate voltage used is the fact that the system 
was let to equilibrate in it. Moreover, it is clear that performing the experiments 'quasi-statically', waiting a long 
enough time at each voltage point, will retrieve the normal field effect. Thus, the voltage scan rate must enter. Fig. 
[7] demonstrates the results of scanning at different rates. Intuitively, relaxing to lower in energy (deeper) metastable 
states should indeed result in a lower conductance. The theoretical framework for understanding this property will 



be given in section ( VI D ) 



The width of the dips (also referred to as cusps) is an important characteristic of the system: an elaborate set of 
experiments |83j have shown that the cusp shape does not depend on the scan rate (for scan rates varying over 3 
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decades), and is independent on disorder and magnetic field. As such, it manifests an intrinsic property of the system, 
which has also been shown to be correlated with the interaction strength (by showing a strong density dependence) |84j . 

We now go on to describe a slightly more involved experimental protocol, that can serve as a useful characterization 
tool to probe changes in the system timescales. 

B. Two-dip experiments 

An important protocol for demonstrating and quantifying the memory effects in electron glasses in a striking way 
was developed by Ovadyahu et al. [HHISS]. The experimental protocol is as follows (see also Fig. [3]): 

1. As in the anomalous field effect, the system is let to equilibrate for a long time (t\ in Fig. [3]), typically of the 
order of day or several days. The gate voltage is fixed at a value V\ during this time. 

2. At time t — 0, the gate voltage is changed to a value Vz- 

3. A time t w later, a scan of conductance vs. gate voltage is made. 

Clearly, if t w = 0, we observe the anomalous field effect mentioned earlier. However, for finite values of the waiting- 
time t w , the system exhibits a striking memory effect: while the dip at the original gate voltage V± is gradually 
erased, a new dip begins to form at the gate voltage V2. For this reason this important protocol is termed the 'two-dip 
experiment'. Fig. [8] shows the result of a typical experiment. 

At some time t w , the two dips will be equal in magnitude. This provides a natural experimental time r characterizing 

Conductance 
measurement 




vg 

FIG. 6: Schematic illustration of the field effect transistor setup. A gate voltage is attached to the sample via a separator, and 
the conductance of the sample is measured while changes in the gate voltage charge the sample, positively or negatively. 
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the system. Experimentally, Ref. [TS] demonstrates that this timescale can be consistently measured in various other 
ways: for example, one may change the value of the gate voltage when the system is close to equilibrium, and measure 
how long it takes for the (logarithmic) change in conductance to decay to half its value at one second. We will discuss 



the physical significance of these timescales in section ( VI D ) , and show why the two approaches give similar results 



In Ref. |52j . the temperature dependence of this timescale was carefully measured. The result appears to be quite 
counter-intuitive: as the temperature increases, and one may expect the dynamics to be faster, r actually increases. 
This was associated with a quantum effect, reminiscent of the anomalous temperature effect which occurs in the 
dynamics of the spin-boson problem for coupling to an Ohmic bath and for a certain regime of parameters |86| . 
However, applying a theory constructed for a two level system to one as complex as electron glass is possibly an 
oversimplification, and this point requires further theoretical attention. Furthermore, one should make sure that the 




FIG. 7: The dependence of the conductance on the gate voltage was checked for various scan rates, in a typical field effect 
transistor setup (see Fig. |6|. In addition to a linear trend, related to the "normal" field effect, there is a surprising symmetric 
component, whose amplitude depends on the scan rate (a). Subtracting the linear trend and normalizing the amplitude, the 
rate dependence is eliminated (b). Figure taken from Z. Ovadyhau, PRB 78, 195120 (2008) with permission from the author. 
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measured timescale is indeed an intrinsic timescale of the system, not associated with t\, as mentioned in section 
pIA l. 

A subtle but important difference between the experiments done on InO and those on granular aluminum |16l 134] , 
relates to the evolution in time of the plateaus outside the dips (cusps). While in InO these are static in time, in 
granular aluminum they show a logarithmic dependence on time as well. Recently, T. Grenet and J. Delahaye showed 
this could be consistently understood in terms of metallic screening in the thin films 87J. 



VI. HOW ELECTRON GLASSES REMEMBER - THEORETICAL PICTURE 
A. Connection between the conductance and the occupation number relaxations 

So far we discussed the relaxations of the occupation numbers. One still has to explain how this influences the 
conductance, which is a much more involved property. Namely, in the analysis we assumed that the relaxation modes 
contribute positively to the conductance. Experimentally, it is indeed clear that the conductance is raised when the 
system is perturbed out of equilibrium. However, it is not obvious theoretically, a-priori, why any perturbation of the 
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FIG. 8: A typical experimental measurement for the two-dip experiment protocol. Initially, the anomalous field effect accounts 
for the dip at zero gate voltage, at which the system was equilibrated. The gate voltage is then shifted to a different value, and 
as the system equilibrates a new dip begins to form at the new value while the initial one is slowly erased. Figure taken from 
A. Vaknin, Z. Ovadyahu and M. Pollak, PRL 81, 669 (1998) with permission from the author. 
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occupancies will tend to raise it. Ref. [SB] shows how this can come about for a non- interacting electronic systems 
(in addition to giving experimental results). In the following, we shall give two explanations of this remarkable 
property, where the interactions play a central role. Before that, we shortly review the properties of the equilibrium 
conductance, which will be important for the discussion which follows. 



1. Miller- Abrahams resistor network 



In a seminal paper, Miller and Abrahams showed that the conductance of a system of non-interacting electrons, 
'hopping' between localized states, can be mapped to the determination of the resistance of a resistor-network, where 
the value of the resistor between site i and site j is given by [50] : 



R ZJ =T/[e 2 ^l (30) 



with 7^ the equilibrium hopping rate from site i to site j, see Eq. (J5J). Detailed balance, assumed to hold at 
equilibrium, ensures us that the resistor is well-defined. The intuition behind this is that well coupled sites will have 
a good conductance between them. As one raises the temperature, the effect on 7^ is exponential, and exceeds the 
linear term in the numerator. For a mesoscopic sample (or for numerical purposes) it might be important to deal 
with the connections to the leads. Ref. |89j shows how to incorporate them into the same framework, with the result 



that the resistor between a site and the leads takes the same form as that of Eq. ( 30 ) . When all the energies involved 



are much larger than the temperature, the rates 7^ can be shown to take the approximate form 



7ij „ « p[ _^ _ |E i -Ml + |E J - T Ml + |E,-E 1 | || (31) 

This is commonly used regardless of the above restriction, erroneously. In Refs. [89j EI] it is shown that by taking 
interactions on the local mean-field level, this formula remains correct, but with Ei renormalized due to the Coulomb 
interactions between sites. 

Determining the resistance of the complete network is still not an easy task. Extensive work has been done on 
this problem using percolation theory methods [91] , some using the concept of percolation in phase-space J26l [27] . A 
direct, 'brute force' method is the numerical determination of the resistance, via Monte-Carlo simulations [53]. In the 
following subsection we show a heuristic approach due to Mott, later extended by Efros and Shklovskii to deal with 
the interacting problem, that finds the leading order behavior. 
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2. Variable Range Hopping 

In 1969, Mott considered the problem of hopping conductance, i.e., the conductance occurring through phonon 
assisted jumps (hops) between localized states. In a brilliant analysis, Mott realized that at low enough temperatures 
it is beneficial for the electrons to hop not to the nearest-neighbor, but to neighbors more further away [22] • This is 
in spite of the exponential-in-distance penalty, in order to reduce the exponential-in-energy penalty. The temperature 
determines the scale for the exponential-in-energy penalty, and thus at extremely low temperature the electron tunnels 



over a large distance, to find a site close to it in energy, see Eq. (31 1 . The result is that: 



a^e-Wr* 1 , (32) 

where d is the dimension of the system, and T ~ with v the density-of-states and £ the localization length 
(taking the dielectric constant to be unity). 



3. Coulomb gap 

About 4 decades ago, it was understood that Coulomb interactions can have profound effects on the density-of- 
states (DOS) of electronic systems near the Fermi energy, namely, leading to the vanishing of the DOS at the Fermi 
energy [33] and the emergence of a 'soft' gap near it [M], suggested by an analysis by Efros and Shklovskii. Since 
then, a large number of analytical [Ml [SHI ED ISS] , numerical [9"6l - U01| and to a lesser extent experimental [1021 1103] 
approaches have dealt with this interesting problem. The current understanding is that there is a power-law DOS at 
equilibrium, with exponent depending on the dimensionality. Since it is a direct result of the Coulomb interactions, 
it is called the Coulomb gap. 

For the current work, we shall not attempt to review the literature on the subject. The only ingredient that we shall 
take away from this for the following, is the existence of a soft-gap near the Fermi energy. Surprisingly, this can also 



be understood from the local mean-field theory discussed in section (IV A I [45 . In two-dimensions, Ref. 89J shows 
that the gap obtained by the local mean-field approach is consistent with that predicted by Efros [28] and Raikh [95] 
by a self-consistent equation approach (taking only single electron transitions into account). Fig. [9] compares this 2D 

2 1 E I 

result with the Efros prediction of -^—i . It would be worthwhile to show how this comes about by analytically solving 
the local mean-field equations, which has not been done to this date. 
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B. Why the conductance increases out-of-equilibrium : Coulomb gap 



It is plausible that when pushed out-of-equilibrium, the Coulomb gap is washed out. It is interesting and experi- 
mentally relevant, to ask how the gap is recovered in time. Experimentally, most of the gap would be recovered in 
very short times, as is analyzed extensively by Tsigankov et al. [104] . At long times, there could be two possibilities: 
the first, is that all metastable states manifest the same Coulomb gap, and that the complete recovery of the Coulomb 
gap to each of them is a slow process, for energies close to the Fermi energy. This approach was supported by Yu 
by generalizing the Efros self-consistent approach to the time-dependent domain [105] . The second possibility is that 
as time goes on, the system probes different metastable states, each of which possseses a slightly different Coulomb 
gap, which 'deepens' as the systems relaxes by going into deeper (lower in energy) metastable states. Work in this 
direction was done in Ref. [104] , 

The idea is now as follows: it is plausible that the DOS near the Fermi energy is crucial for the conductance 
properties. This was used by Efros and Shklovskii in a seminal work , to show how the Mott picture for variable- 



range-hopping, briefly discussed in section (VI A 2), should be modified [94. Although the use of the single particle 



tunneling DOS in the conductance calculation is not justified [106 , their prediction was confirmed in a huge number 
of experimental systems [107] , functionally but not necessarily quantitatively. In a recent work [89 , it was shown 




FIG. 9: Comparison between the local mean-field approach [3^ and the Efros self-consistent equation approach 
a linear Coulomb gap (solid line) in two-dimensions, with slope ~ -^j . The figure is taken from Ref. [89] . 



yielding 
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that by following the local mean-field approach, one can self-consistently take the interactions into account, and 
essentially verify the Efros-Shklovskii picture. Thus, the Coulomb gap reduces the conductance, by reducing the 
density of available states near the Fermi energy. The picture that emerges is that through the Coulomb interactions 
the system 'digs' a hole in the DOS, and makes its own conductance diminish. Therefore, by kicking the system out 
of equilibrium, we enhance the conductance, as we assumed earlier in the analysis of the IRM protocol, by taking the 
contribution of the modes to be positive. The problem that remains is the understanding of the timescales: as we 
mentioned, various references suggest that the Coulomb gap is created on fast timescales of the order of the Maxwell 
time [91] (although, as mentioned, some works indicate that at energies close to the Fermi energy, i. e. , the bottom of 
the gap, the Coulomb gap takes a long time to form [105J ) . The next section will give a possible solution, that shows 
how this effect can be understood on more general ground, beyond the single particle relaxations. 

C. Why the conductance increases out-of-equilibrium : long time relaxations 

In the presence of Coulomb interactions, slow relaxing modes can modify the energies of other sites, and thus 
influence the conductance even if the modes are isolated, and the tunneling between them and the rest of the sites 
is negligible (since it is exponentially suppressed). This picture also explains why one may talk about an out-of- 
equilibrium conductance: there can be a separation of timescales between the (relatively quick) equilibration time of 
the current carrying sites, the 'backbone', and the isolated, slowly relaxing modes outside it. When the slow modes 
have not equilibrated with the rest of the system, strictly speaking the system is out-of-equilibrium. Nevertheless, 
there is a well defined conductance that can be experimentally measured, defined by the internally equilibrated sites 
belonging to the backbone, which have reached quasi-equilibrium in the presence of the slowly changing potential 
created by the slow modes (i.e., their statistics is described by the Fermi-Dirac distribution). In this sense, the 
combination of the long ranged Coulomb interactions and the localized modes explains why the conductance returns 
slowly to its equilibrium value. This is not enough, however, to explain why it does so monotonically, as is clearly 
seen experimentally in the IRM experiment, for example. 

In order to understand this, a simple argument has been put forward by Ref. ;29J. In the following we give our 
version of the argument. 

An essential ingredient is the separation of timescales between those associated with hopping between sites which 
contribute significantly to the conduction (by definition "fast" processes), and clusters of sites which support "slow" 
modes: Ref. |20| shows that isolated clusters lead to an abundance of slow modes with a distribution approximately 
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given by P(A) ~ 1/A, but this is not the only mechanism leading to slow modes: Coulomb interactions will give rise 
to modes related to the tunneling of many electrons, related to the work of Kozub et al. [35]. The contribution of 
these clusters to the conductance is not negligible, however, since they affect the energies of the sites associated with 
the fast processes (from now on we refer to these as the "backbone"), through the Coulomb interactions. One might 
expect that the contribution will be effectively random: sometimes they will give rise to an increase and sometimes 
a decrease. It turns out that on average, however, they give rise to a positive contribution. We now go on to the 
crux of the matter: let us consider a relaxation mode associated with an isolated cluster some distance away from the 
backbone, as shown in Fig. flO| 

Since we have a separation of timescales, it is still justified to think of the Miller-Abrahams resistance network 



associated with the conduction backbone (see section (VIA1), which is presumably in equilibrium, but with the 
energies of the sites in it affected by the current configurations of the slowly relaxing modes outside the backbone, 
which have not yet fully equilibrated. Let us consider the value of the resistance associated with two close sites on 



the backbone, A and B in Fig. 10 The resistance between them is given by Eq. (31 1, showing that the dependence 
on energy is: 
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FIG. 10: Schematic illustration used to demonstrate why a relaxation of a mode outside the conduction "backbone" contributes 
negatively to the conduction, on average. The dotted (purple) arrows denote the dipole moments of the cluster and of two sites 
in the backbone. These will prefer to point in opposite directions, on the average, and thus while a cluster relaxes it will tend 
to raise the energy difference between the two sites in the backbone, on average, thus making the conductance lower. 
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r> B JEa-mI + IEb-^I + IEa-Eb^ f . 

R ~ exp[ — j (33) 

The claim is that modes that relax in energy will tend to make the value of the resistor larger (on average) as they 
relax. Let us assume that Ej± < fJ- and Eg > fi. We know that in equilibrium (which is assumed for the sites on 
the backbone), the average occupations obey Fermi-Dirac statistics [15], therefore we know that site A is negatively 
charged and site B is positively charged. Therefore there is a net dipole moment with a component in the positive 
x direction, see Fig. |10[ On average, it is plausible that there will therefore be a tendency for the dipole moment of 
the relaxation cluster (encircled by an eclipse in the figure) to be aligned in the —x direction. Therefore, on average, 
in the relaxation process the cluster will contribute to the electric field in the vicinity of the sites A and B in the 
positive x direction, thus making the energy difference between them larger still, and increasing the resistance of Eq. 



(33). It is possible that the argument does not rely explicitly on the specific form of the interactions, as long as it 
is not too short ranged. Thus, it seems like a robust and generic mechanism (albeit, non-rigorous in the form given 
here) for explaining why the conductance decreases as the system relaxes, which can be used to explain qualitatively 
some of the experiments described earlier. 

D. Anomalous field effect and two-dip experiments - qualitative theory 

Following the theoretical discussion, we can go back to explain qualitatively the physics behind the anomalous field 
effect and the two-dip experiment. 



It is clear from the argument of section (VIC), that the conductance decreases as the system equilibrates. Already 



after single particle transitions take place, the Coulomb gap is formed, as discussed in section (VIB) which accounts 
for most of the reduction in conductance. Nevertheless, as many particle rearrangements take place, the system 
conductances continues to decrease, which explains the long timescales involved. 

What happens when we make a gate voltage scan? Let us consider a "fast" scan, without yet mentioning compared 
to what. It is clear that at a given point of the scan, the system is nearly at equilibrium with respect to the gate voltage 
at which it was equilibrated, not at the current value of the gate voltage. Thus, the system is out-of-equilibrium with 
regard to the current gate voltage, and therefore has an improved conductance. For small changes of gate voltage, 
such that the system configuration is still approximately in a local minimum (metastable state), this statement is not 
true: for these, the system is nearly in equilibrium, and the value of the conductance is nearly the lowest one possible. 
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This explains qualitatively why one obtains an anomalous field effect, without explaining the characteristic scales 
involved, the dependence of the scan rate, and without considering the normal field effect which will give an additive 
linear dependence to the anomalous effect. 

It should be emphasized that the Coulomb gap gives intuition of how interactions diminish the conductance of the 
system, but looking at the time scales involved it is more plausible that it is formed "immediately" at each point of 
the scan, and that the anomalous field effect is due to more than single particle transitions (as mentioned earlier, Ref. 
shows that the Coulomb gap due to single electron tunneling forms in very short times). 



Let us proceed to the explanation of the two-dip experiment, described in section ( V B I, which is relatively straight- 
forward on the qualitative level once the anomalous field effect is understood. If we make a gate voltage scan immedi- 
ately after the gate voltage is changed to its new value, clearly we are still measuring the anomalous field effect, since 
the system did not change its configuration yet. The result of a scan made after a very long time is also clear: the 
system will equilibrate at the new gate voltage, and we will measure an anomalous field effect shifted to the new value 
of gate voltage. Thus, over time the initial dip must slowly vanish, and a new one must begin to form. It turns out 
experimentally that in many cases there is a symmetry between the depth of the new dip formed and the amplitude 



by which the old dip has become smaller. This can be readily understood based on the results of section (IVA2|: 
in a two-dip protocol, the system is initially in equilibrium, and then the gate voltage is changed to a new value. 
Let us assume a quick scan is made after time t w . If we take from this scan only one value, the conductance at the 
new gate voltage V2, which is exactly the position of one of the two dips, we obtain a similar dependence to that of 
the IRM protocol on time during stage II of the IRM protocol (the only difference is that here there is an additional 
perturbation made, during the time of the scan itself, t scan , which slightly destroys the new dip formed, and has to 
be accounted for |108j V We therefore find that the new dip is formed with an amplitude growing logarithmically in 
time. Let us denote by tdi P the time it takes to scan the width of the cusp. Therefore measuring the conductance at 
the old gate voltage in the two-dip experiment is approximately equivalent to the IRM protocol in Stage III, with a 
relevant time of the order of tdi P - this is the time the system effectively had to dig the gap at V\. Thus, the result 
should be log(l + t w /td ip ), and for t w 3> t^ip the formation of the new dip and the erasure of the old one both scale 
as log(t w ), which accounts for the symmetry described above. Now we are in a position to understand the timescale 



involved in the two-dip experiment, mentioned in section (VB|: the requirement that the dips are of the same depth 
reduces to the equation: 
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\\og(r /t d ip)\ = log(r/*di p ) = |log(A min r) - \og(t scan /t dlp )\ = -\og(X min Tt dip /t scan ), 

where the term — log(t scan /tdi P ) accounts for the reduction of the new dip formed due to the scan. 
This finally leads to the relation: 



(34) 



Amir 



(35) 



Typically, t scan is of the order of tens of seconds [53]. Thus measuring a value of r of the order of 1000 seconds 
implies an associated timescale 1/A m ; n of about a day. This is possibly due to the finite amount of time the system 
was let to equilibrate prior to the experiment, and not the true cutoff of the underlying distribution. 

Before continuing with the two-dip experiment, it is useful to discuss now another way of measuring this timescale, 
which experimentally leads to similar results [15J . and we shall show that it is also consistent theoretically: in this 
method, one changes the gate voltage, and then simply waits until the conductance deviation decays to its half value 
from that measured after one second. 

The procedure is demonstrated in Fig. |11| 



Here, the analysis is straightforward: according to Eq. (14) the decay in this procedure goes as —je ~ log(A m i n i), 
thus, the normalized decay is: 
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FIG. 11: Illustration of one of the procedures that can be used to determine the timescale of the system: the system is let to 
equilibrate for a long time, after which a sudden change in gate voltage is made. The signal (conductance change) is normalized 
to unity at t = 1 second, and the time at which the (logarithmically relaxing) conductance reaches the value of half defines the 
timescale of the system. Eq. (1351) shows the physical significance of this timescale. 
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Sa(t = lsec) -j E - log(A min ) 
Equating this to 1/2, we find that: 



5v(t) _ -Je - log(A min t) 



^min 

Clearly, due to the choice of normalization after one second, similar to the value t s , the two methods qualitatively 
agree. 



Let us return to the two-dip experiment. The theory described in section (IV A3 1 might hint that one can use 
linear response to understand also the two-dip protocol: if the response to a generic period sequence of square pulses 
can be decomposed into a superposition of contributions arising from the changes in gate voltage, i.e., 6a(t) ~ 
/ G(t — t')^r(t'), where G{t — t') is a function describing the response of the system to a step function in gate voltage, 
why not plug in a more complicated dependence of gate voltage on time, V(t), to understand two-dip experiments 
quantitatively? This approach is incorrect, due to a fundamental reason: the dependence on gate voltage is not linear. 
For square pulses, only two values of gate voltage play a role, and therefore the assumption of linearity is unnecessary, 



which is why the analysis of section (IV A3) is correct. But for the two-dip experiment, if we are interested in the 
values of the conductance associated with other values of gate voltage than V\ or V2, we cannot use this linear response 
equation. 

Therefore within the framework of the 1/A relaxation rate distribution we cannot explain quantitatively the shape 
of the dips, as well as other interesting questions regarding it, e.g: What determines their dependence temperature 
and disorder? 

Various works have addressed these questions, which we briefly mention now. The general physical picture presented 
here is similar to that of Rcfs. 29, 109J with regard to the mechanism leading to the decrease of the conductance 
and the anomalous field effect, which is the basis for understanding (qualitatively) the two-dip experiment. Ref. [110] 
suggests that the slow relaxing modes are not (as was assumed above) an intrinsic part of the electronic system, but 
rather, atomic configurations which change over time and interact with the electronic configuration. Ref. [29] , by the 
same authors, explain that the same mechanism would also be consistent with an intrinsic picture, similar to that 
described here, and refer to a specific subset of relaxation modes, with a definite structure in space, which they refer 
to as "chessboard cluster". These were introduced by the same authors in Ref. |111) in order to explain the 1// 
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noise which occurs in electronic glasses, which we discuss in section (IVB 1 ). Ref. [112] also deals with the two-dip 
experiment, but focuses more on the dependence of the dip on temperature, carrier density and disorder, and less on 
the timescales associated with the formation of the phenomenon. 



E. Estimating the timescales 

So far we have mainly discussed the functional form of the relaxations involved, but have not demonstrated why 
the long time cutoffs of the underlying distribution exceed the experimental timescales involved (which can be hours 
and days). In fact, this remains a difficult question, for which no clear cut answer exists. While many theoretical 
mechanisms that can give rise to such slow relaxations exist, it is not clear what distinguishes the slow relaxations 
described in the previous section, on Anderson insulators and granular systems, and various experiments performed on 
semiconductors, for which the experimental results are different, and no slow relaxations with full aging is observed. 
An experimental clue in the direction of solving this profound problem is given by Ref. [84] . suggesting that the 
electron density plays a role in determining the system relaxation time. Indeed, the electronic density in doped 
semiconductors below the metal-insulator transition does not reach the values of indium oxide. Recent experiments 
on semiconductors shows that these systems may also exhibit rich behavior and non-exponential relaxations, albeit 
on much shorter timescales [1 13] . Measuring the density dependence of the relaxation timescales involved might shed 
new light on this problem. For the indium oxide samples, on the other hand, it would be good to check experimentally 
if for the low-density samples one can see a deviation from the clean logarithmic relaxations, signalling that we are 
nearing the cutoff of the relaxation rate distribution. An idea suggested by Ovadyhau [52] emphasized the effect of 
Anderson orthogonality catastrophe is slowing down the rate: indeed, experiments suggest that the slow relaxation 
regime corresponds to that where more than a single electron is present in a localization volume [3TJ II 14] . Another 
ingredient which may be important for understanding the source of the slow relaxations is the simultaneous tunneling 
of more than one electron [21H291 11151 1116] . However, such many-particle tunneling events can be present also in 
semiconductors: one has to show that their significance depends on a parameter which is different in the two systems 
- a natural candidate is the ratio of the localization length to the average nearest-neighbor distance. This criterion is 
equivalent to the above mention criterion of having many particles in a localization. Showing this systematically by 
analyzing the statistics of the many-particle tunneling processes would be a significant step. 
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VII. FUTURE PROSPECTS 

In this review we discussed the current experimental and theoretical understanding of the memory and aging 
associated with electron glasses. Throughout, a number of open questions, mostly theoretical, were mentioned. 
Perhaps the most fundamental regards the timescales associated with the relaxations. What is the role of many 
particle transitions in determining it? How many electrons tunnel simultaneously? Does the Anderson Orthogonality 
Catastrophe (or a FranckCondon type effect) play a substantial role? 

Another deep question regards the basic ingredients needed to observe the form of aging and slow relaxations 
discussed in this review. Experiments performed on spin glasses, for example, show more complex aging behavior 
II IT] . It is in our opinion a worthwhile question to understand the similarities and differences between these two 
related systems. Structural glasses also show similar behavior |118j . which should be further explored. 

More technical questions, to which at present we do not have a full answer, regard the quantitative understanding 



of the two-dip experiment (section (VB|. Also, the TRM protocol described in section ( III B I still has to be explored 
further. 

Altogether, it seems that electron glasses still pose many more questions and challenges, and prove a useful platform 
for investigating glassy behavior. 
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